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Dark Periods in Rabi Oscillations of Superconducting Phase Qubit Coupled to a 

Microscopic Two-Level System 
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We propose a scheme to demonstrate macroscopic quantum jumps in a superconducting phase 
qubit coupled to a microscopic two-level system in the Josephson tunnel junction. Irradiated with 
suitable microwaves, the Rabi oscillations of the qubit exhibit signatures of quantum jumps: a ran- 
dom telegraph signal with long intervals of intense macroscopic quantum tunneling events (bright 
periods) interrupted by the complete absence of tunneling events (dark periods). An analytical 
model is developed to describe the width of the dark periods quantitatively. The numerical sim- 
ulations indicate that our analytical model can capture underlying physics of the system. Besides 
calibrating the quality of the microscopic two-level system, our results have significance in quantum 
information process since dark periods in Rabi oscillations are also responsible for errors in quantum 
computing with superconducting qubits. 

PACS numbers: 74.50.-|-r, 85.25.Cp 
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I. INTRODUCTION 



Superconducting Josephson devices coherently driven 
by external fields provide new insights into fundamen- 
tals of quantum mechanics and hold promise for use in 
quantum computation and quantum informationfii^i^ Re- 
cent experiments based on Josephson tunnel junctions 
have unambiguously demonstrated the quantum behav- 
ior of macroscopic variable i^i^ii One interesting quan- 
tum phenomenon is known as quantum jumps, which 
was proposed by Bohr as early as 1913.^ Bohr suggested 
that the interaction of light and matter occurs in such 
a way that an atom undergoes instantaneous transitions 
of its internal state upon the emission or absorption of 
a light quantum. These sudden transitions have become 
known as 'quantum jumps'.— 'i^'^'i^'^ Quantum jumps 
were firstly observed in experiments in the 1980s by in- 
vestigating the fluorescence of a single trapped ion driven 
by laser In such atomic systems, quantum jumps are a 
random telegraphic process with long intervals of intense 
photon emissions interrupted by periods of the absence of 
photons. To observe quantum jumps in macroscopic sys- 
tems such as superconducting qubits, one may naturally 
think that it can be achieved by detecting microwave pho- 
ton emissions in analogy with the way in atomic systems. 
However, this idea suffers from the absence of microwave 
photodetectors, although some theoretical efforts have 
been madeJ^ii^ Recently, macroscopic quantum jumps 
were experimentally demonstrated for the first time in a 
superconducting phase qubit coupled with a TLS inside 
the Josephson junction. In this experiment, the state of 
the system is read out not by detecting photon emissions 
but by detecting macroscopic quantum tunneling events, 
and quantum jumps behave in the form of jumping ran- 
domly between upper branch and lower branch of the 
switching currents. In addition, recent experiments on 
superconducting charge qubits also demonstrated quan- 



tum jumps by observing quasiparticle tunneling in the 
time domain ) ^^d^ and theory for the kinetics of the sys- 
tem has been developed^ 

In this article, we show that macroscopic quantum 
jumps can be better observed in a Rabi-oscillation exper- 
iment in the phase qubit-TLS coupling system. The main 
points we are going to make are the following. Firstly, 
quantum jumps in coherent excitations of macroscopic 
quantum states can be well studied based on our model. 
In the experiment by Yu et al.^ the energy level struc- 
ture keeps changing during the measurement time, and 
coherent characteristic of the dynamics cannot be ob- 
served directly. The scheme proposed in this article is 
based on a fixed energy level structure, which is usually 
implemented to demonstrate Rabi oscillations as the ba- 
sic skill to manipulate quantum states. As we can see 
below, as a result of the coherent dynamics in our model, 
some new features of quantum jumps which have not 
been expected in the previous experiments can appear. 

Secondly, macroscopic quantum jumps in the qubit- 
TLS coupling system bring new insights into our un- 
derstanding of the effects of TLS on the dynamics of 
the qubit. Qubit-TLS coupling system has received dra- 
matic attention recently in both experimenta P^'^^i^^ and 
theoretical— i^^'^^i^'^'^^i^^ studies because it is suggested 
that TLS is a major source of decoherence in supercon- 
ducting Josephson qubits. In particular, Rabi oscilla- 
tions in such qubit-TLS coupling system have been stud- 
ied by several authors i^ii^liS^iS^ However, the previous 
works mainly emphasized on the ensemble characteristic 
of the system, without considering the effects of quantum 
jumps in single trajectories, which is actually the origi- 
nal form of experimental data.^'^^ Therefore, a detailed 
description of the effect of TLS on the trajectories of a 
single qubit is desirable. 

Thirdly, recent experiments have demonstrated that 
TLS can be used as quantum memory with good 
performance.'^" It is also suggested that such TLSs them- 
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selves can serve as qubits^i^ and can be implemented to 
generate genuine multi-qubit entangled statesi^ There- 
fore, characterizing the TLS is critical to improve the 
performance of such solid-state qubits. In this article, we 
show that the macroscopic quantum jumps phenomenon 
can be used to read out the state of TLS, which provides a 
new approach to calibrate individual TLS quantitatively. 

This article is organized as follows. In Sec. II we briefly 
describe the basic physics of superconducting phase qubit 
under coherent driving fields, and then introduce the 
Monte Carlo wavefunction method which can be adopted 
to study the Rabi oscillations in the superconducting 
phase qubit. In Sec. Ill we introduce the qubit-TLS cou- 
pling system, and present the Hamiltonian of the hybrid 
system. In Sec. IV, we provide a clear physical picture of 
macroscopic quantum jumps in such qubit-TLS coupling 
system. In Sec. V, we study the width of dark periods in 
Rabi oscillations quantitatively with both numerical and 
analytical methods. In Sec VI, we show that the quan- 
tum jumps approach can be used to characterize the TLS 
quantitatively, and this article ends with a brief discus- 
sion on the relation between quantum jumps and readout 
fidelity in superconducting quantum computing. 

II. QUANTUM JUMP APPROACH TO RABI 
OSCILLATIONS IN SUPERCONDUCTING 
PHASE QUBIT 

In this section we briefiy present the basic physics of 
superconducting Josephson phase qubit driven by coher- 
ent fields, and then introduce how the quantum-jump 
approach can be used to simulate Rabi oscillations of a 
single qubit. To grasp the spirit of quantum-jump ap- 
proach in a simple way, here we firstly consider the qubit 
without coupling to a TLS. 

Superconducting Josephson phase qubit is essentially 
a current-biased Josephson junction. The Hamiltonian 
of the phase qubit as shown in Fig. 1(a) reads^i^i^i^l 

Hqb^ T^Q - ^—cosS- —6, (1) 
zu Zn Zn 

where Iq is the critical current of the Josephson junc- 
tion, / is the bias current, C is the junction capacitance, 
$0 = h/2e is the flux quantum, Q denotes the charge 
operator and 5 represents the gauge invariant phase dif- 
ference across the junction, which obeys the convectional 
quantum commutation relation [S, Q] = 2ei. Quantum 
behavior can be observed for large area junctions when 
the bias current is slightly smaller than the critical cur- 
rent. The junction works as a phase qubit when the 
Josephson coupling energy Ej — /o<I>o/2vr is much larger 
than the single charging energy Ec — e^/2C. In this 
regime, the two lowest energy levels, |0) and |1) as shown 
in Fig. 1(b), are usually employed as a qubit in quantum 
computation. The state of the qubit can be controlled 
through the bias current I{t) given by^^ 

I{t) = he + ^I{t) = he + hw cos Lut, (2) 
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FIG. 1: (a) Equivalent circuit of a current-biased Josephson 
tunnel junction, (b) An illustration of various coherent and 
incoherent processes in a superconducting phase qubit radi- 
ated by a microwave. 

where the classical bias current is parameterized by the 
dc current he and the ac current with the magnitude /^^ 
and frequency to. 

Truncating the full Hilbert space of the junction to the 
qubit subspace {|0), |1)}, the Hamiltonian of the phase 
qubit can be written as 

Hgb^hcuo\Q){Q\ + 'huji\l){l\ + nn^ cosc^t(|0)(l| + |l)(0|), 

(3) 

where w„ is the energy frequency of state |n), Qm = 
hw '\/l/2?ia;ioC is Rabi frequency, and wio = toi — loq is 
the energy frequency between state |0) and |1). In the 
interaction picture and choosing a rotating frame of the 
frequency w, Hamiltonian ^ can be simplified to 

jy,, = ?iA|l)(l| + ^(|0)(l| + |l)(0|), (4) 

where A = ujiq — lo represents the detuning. 

To better understand our simulation method, we give 
a brief description of the experimental procedure firstly. 
In experiments;^ the system is prepared in the initial 
state |0) firstly, then a microwave source is turned on 
for a duration time (For convenience, we call the 
measurement time below). Because the tunneling rates 
depend exponentially on the barrier height, the bias cur- 
rent Ida can be chosen so that the tunneling from |0) is 
essentially 'frozen out':^ Therefore, the Rabi oscillations 
between states |0) and |1) lead to an oscillating proba- 
bility for the system to tunnel out of the potential. At 
the end of the measurement time for a single run, no 
matter the tunneling event has happened or not, the bi- 
ased current he is adjusted to initialize the state for the 
next run. Here we label the time interval for initializing 
the qubit state as Tp. 

To simulate the dynamics of the phase qubit system 
with and without TLS coupling, we adopt Monte Carlo 
wavefunction method, which is also called 'quantum tra- 
jectories' method^iiSiiiii^i^i^ada. In this method, the 
time evolution of the system can be described by the 
non-Hermitian effective Hamiltonian 

i/e//=i?,6-y(7io + ri)|l)(l|, (5) 
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FIG. 2; (Color online) Simulated Rabi oscillations of a single 
superconducting phase qubit. (a) Average over A'^ = 5000 
quantum tunneling events obtained using Monte Carlo wave- 
function metliod(green dots). The red line is an ensemble 
result obtained using master equation, (b) Trajectories of 
quantum tunneling events obtained with Monte Carlo wave- 
function method. The parameters are Qrn = lOFi, A = 0, 
and 710 = Ti/A. 

where 710 is the rate of energy relaxation from |1) to |0) 
and Fi is the tunneling rate from the level |1) out of the 
potential. The 'conditional state' |^'c(0) of the system 
is determined by 

\^cit + St)) =CXp[-lHeffSt/h]\^,{t)). (6) 

The procedure adopted in quantum-jump simulation of 
Rabi oscillations of a single qubit can be summarized as 
follows: 

(i) Determine the current probability of a relaxation or 
tunneling event, i.e., APr = jioSt{'^\l){l\^) and APf = 

ri6t{^\i){i\^). 

(ii) Obtain random numbers ri and r2 between zero 
and one, compare with APt and APr respectively, and 
decide on tunneling or relaxation for different cases: 

(a) If Ti < APf a quantum tunneling event happens. 
Register the time t of the event, and then turn to step 
(in). 

(b) If ri > APt and r2 < AP( there is a relaxation, so 
that the system jumps to the state |0). 

(c) If ri > APr and r2 > APt no quantum jumps take 
place, so the system evolves under the influence of the 



non-Hermitian form 

{i-m)H,t5t^{j/2)st\i){i\m 

' (l-AP)i/2 ' ^'^ 

where 7 = 710 + Fi and AP = AP^ -I- APt. 

(iii) If no tunneling event happens, repeat the previous 
steps until the end of the measurement time Tm- 

(iv) Accounting time t over many simulation runs. 
By using this Monte Carlo wavefunction method, we 

simulate the trajectories of Rabi oscillations in a single 
superconducting phase qubit and the results are shown 
in Fig. 2(b). We emphasize that it is possible that a tun- 
neling event never happens during a finite measurement 
time Tm for a single run, which is important for observing 
dark periods as discussed in Sec. V. To test the efficiency 
of our simulation method, we compare the simulated re- 
sults with those obtained with master equation used in 
Ref [6]. As shown in Fig. 2(a), they agreed very well. Ac- 
tually, it can be proved that the Monte Carlo method 
is equivalent, on average, to the master equatioui ^i^^i^^ 
Guaranteed with the efficient method for simulating Rabi 
oscillations in superconducting phase qubit, we turn to 
the more interesting qubit-TLS coupling system. 

III. HAMILTONIAN OF QUBIT-TLS 
COUPLING SYSTEM 

Experiments have shown that some TLSs may locate 
inside the Josephson tunnel barrier, as illustrated in 
Fig. 3(a). A TLS is understood to be an atom, or a 
small group of atoms, that tunnels between two lattice 
configurations f^Siii^ with different wave functions \L) and 
\R). The two states of the TLS correspond to two dif- 
ferent values of the Josephson junction critical current 
Iq which is proportional to the square of the tunneling 
matrix element across the junction. When the TLS is in 
state \R) (|P)), the junction critical current is Iqr {Iql)- 
Then the interaction Hamiltonian between the qubit and 
the TLS is: 

H,nt = -^^cos<5®|i?)(P|-%^cos5®|L)(L|. (8) 

Assume an asymmetric potential with energy sepa- 
rated by Tildtls for the TLS, then the ground and ex- 
cited states are \g) = sin(6'/2)|L) +cos(0/2)|P) and |e) = 
cos(6l/2)|L) - sin(6l/2)|P), where uotls = V^^ + A^, 
~ arctan(Ao/e) with e/2 being the energy asymme- 
try and Ao/2 being the bare tunneling matrix element 
(Fig. 3(b)). Considering the junction is biased near its 
critical current^^i^i, i.e., 5 — tt/2 — 5' with \5'\ <C 1, 
then cos(7r/2 — 5') = sin 5' w 5' . In this case, 5' can be 
well approximated as position coordinate operator of the 
harmonic oscillator. Using these facts and including only 
the dominant resonant terms arising from the interaction 
Hamiltonian, Equation ([5]) becomes 

H^r^t = ^^^y^^(|0e)(l5| + |lff)(Oe|), (9) 
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FIG. 3: (Color online) (a) Schematic of a TLS in the oxide 
tunnel barrier of the Josephson junction. Because the barrier 
material is disordered, some atoms can occupy two positions, 
labelled \L) and |7i). (b) The position states are separated by 
an energy difference e and connected by a tunneling energy 
Ao/2. In the energy eigenstate basis, the ground state \g) 
and excited state |e), are separated by an energy ujtls- The 
charge motion between \L) and |_R) couples to the currents 
and voltages of the qubit circuit. 
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FIG. 4: (a) V^-type energy structure for observing quantum 
jumps in the fluorescence of a single atom. Two upper levels 
jl) and 1 2) couple to a common ground state |0). The |1) ^ 
|0) transition is assumed to be strong while the |2) ^ |0) 
transition is weak, (b) Schematic energy level diagram for 
a junction coupled to a TLS. The ground state (the excited 
state) of the qubit is denoted as |0} (|1}); \g) and |e) represent 
the ground state and the excited state of TLS. The \lg) ^ 
|0e) transition is weak, while the |0e) <-> |le)(|0(?) ^ llff)) 
transition is strong. 



where SIq = Iqr — Iql is the fluctuation amplitude in 
Iq produced by the TLS. The coupling of the two in- 
termediate energy levels through Hint produces an en- 
ergy sphtting which can be characterized in spectroscopic 
measurementsJi2ii^i22, The Hamiltonian of the qubit- 
TLS couphng system in the basis {\0g), \^g), |0e), and 
|le)} then reads 



H = h 



flm cos Lot 

r2„i COS ujt ujiQ ric 

ilc I^TLS ilmCOSUjt 

iljri COS Ujt U) 10 + ujtls 



(10) 

To understand the underlying physics that governs the 
dynamic of the system more clearly, we chose the inter- 
action picture and make a transformation to a rotating 
frame. Then the Hamiltonian can be simpHfied to the 
time independent form (see appendix) 



H' 



r2™/2 

ric A + Ar ^„i/2 

n„j2 2A + Ar 



(11) 



where A = wio — lo and A^ = ojtls ~ ^lo are the de- 
tunings, Q.C — (SIq sm0/2)y^l/{2fiujiQC) i s the avoided 
energy level crossing between |0e) and \lg), and fJ,,,. = 
Ifj^w \/^/2huJiQC is the Rabi frequency between the qubit 
state |0) and |1). From Hamiltonian (fTT|) . we can di- 
vide the system into two subspaces A = {\0g), |lg)} and 
B = {|0e), |le)}. The effective coupling strength between 
A and B is decided by several parameters including Qc, 
Ar, A and O^. Interestingly, the parameters A^, A and 
flm can be easily controlled by adjusting the bias current 
Idc, the microwave frequency and the microwave ampli- 
tude, respectively. 1^ While in the interior of the subsys- 
tem, the coupling strength ^lm/2 between states \0g) and 



1 1.9) (|0e) and |le) ) can be controlled by adjusting the 
amplitude of the microwave.^^ Therefore, the qubit-TLS 
coupling system behaves as an "artificial atom" which 
can be manipulated flexibly. 



IV. QUANTUM JUMPS IN QUBIT-TLS 
COUPLING SYSTEM 



To have a clear physical picture of macroscopic quan- 
tum jumps in qubit-TLS coupling system, we provide an 
analogy with quantum jumps in quantum optics, which 
has been well developed for years -^iio^iiii^ii^iiii^SiSQ^ 

quantum optics, the system generally has energy struc- 
ture shown in Fig. 4(a). Two excited states |1) and |2) 
are connected to a common ground state |0) via a strong 
and weak transition, respectively. The fluorescent pho- 
tons from the strong transition are observed. However, 
an excitation of the weak transition where the electron is 
temporarily shelved in the metastable level |2) will cause 
the strong transition to be turned off. Therefore, it is 
possible to monitor the quantum jumps of the weak tran- 
sition via the signal provided by the fluorescence of the 
strong transition. In the language of quantum measure- 
ment theory,*^ the fluorescence from the strong transition 
acts as a pointer from which the microscopic quantum 
state of the atom may be determined. 

Similarly, in the qubit-TLS coupling system consid- 
ered here, quantum jumps between macroscopic quantum 
states are also proposed to happen. As discussed in Sec. 
HI, we can have a weak transition between states \lg) 
and |0e) by adjusting the bias current and a strong tran- 
sition between |0e) and |le)(or between \0g) and \lg)) 
by adjusting the microwave amplitude. However, differ- 
ent from observing fluorescent photons in atom systems, 
the states of the qubit-TLS coupling system are read 
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FIG. 5: (Color online) Dark periods in Rabi oscillations of 
qubit-TLS coupling system. The detunings are (a) = 2Tie 
and (b) Ar — SFie. Lower panel: the filtered data corre- 
sponds to on and off of the Rabi oscillations. The param- 
eters for simulation are A = 0, 710 — rie/4, flc = rie/4, 
fim = lOFie, and = 10/Fie. 



out through quantum tunneling. As has been proved 
in experiments^'', in the qubit-TLS coupling system, dif- 
ferent states correspond to different tunnehng rates .~ 
Therefore, by adjusting the bias current at a proper 
value, only the tunneling from state |le) is prominent, 
and tunnehng from other states can be neglected. 

Then the physics for quantum jumps in qubit-TLS 
coupling system is straightforward. When the sys- 
tem stays in subspace i3{|0e), |le)}, quantum tunneling 
events in the Rabi oscillation form (see Fig. 2) can be ob- 
served. However, once the system transitions into sub- 
space ^{lOg), (dark area in Fig. 4(b)), where the 
system has little probability to tunnel out the potential, 
no quantum tunneling events can be observed during the 
measurement time, and then dark periods in Rabi oscil- 
lations appear. To substantiate our proposal, we investi- 
gate this phenomena numerically with the Monte Carlo 
wavefunction method introduced in Sec. II. In the qubit- 
TLS coupling system, the dynamics can be described by 
the non-Hermitian effective Hamiltonian 



H, 



^'-?(7io + ri.)|le)(le|-^7io|l5>(l3l, (12) 



where Fie is the tunneling rate of state |le) and 710 is 
the relaxation rate from |le) to |0e) (|l.g) to \0g)). It is 
noticed that we do not take into account the relaxation 
effect of the TLS, because the lifetime of the TLS is much 
longer than that of the phase qubitj ^^'^^i'^^ Actually, this 
effect can be easily considered in the same way as we 
consider the relaxation of the qubit. The quantum state 
either evolves according to the schrodinger equation, or 
'jumps' to an eigenstate of the system with certain prob- 
ability. The measurement result serves as a pointer to 
determine the quantum state of the system. Once a quan- 
tum tunneling event happens, we know that the state of 



FIG. 6: (Color online) The time evolution of the population of 
state I le) . The dark periods correspond to the shelving of the 
system in subspace A{|0(;), The parameters are A = 0, 

710 = ri./4, = Fie/4, = lOFie, and A, = 2Fie 



TLS is |e); then in the next run the initial state of the sys- 
tem is prepared in the corresponding qubit ground state 
|0e). On the contrary, if no quantum tunneling event hap- 
pens during a sufhcient long measurement time Tm, we 
can confirm that the system is in subspace yl{|05), \ lg)}, 
i.e., the state of TLS is \g). Then in the next run of sim- 
ulation the initial state is prepared in the corresponding 
qubit ground state \0g) and dark periods in Rabi oscilla- 
tions may appear. We emphasize that if the measurement 
time T„i is short, the TLS's state cannot be determined 
if no quantum tunneling event happens at the end of Tm- 
In this case, we decide the initial state of the next run 
using Monte Carlo method according to the population 
of \g) and |e) respectively. 

As shown in Fig. 5, in the simulated trajectories the 
quantum tunneling events are disturbed by dark periods 
during which no tunneling events are observed. In ad- 
dition, it is found that the average width (Nd) of dark 
periods for — 2Fie is smaller than that for = SFie 
(Here (No) means the the average number of runs in a 
single dark period). One may simply think that as the 
detuning A^ increases, the coupling between subspace A 
and subspace B becomes weaker. Therefore the system 
is more difficult to make transitions between A and B, 
thus leading to a larger average width of dark periods. 
However, this is not always the case. As discussed below, 
we will give a quantitative description between (No) and 

Ar. 



V. RELATION BETWEEN DARK PERIODS 
AND DETUNING 

In this section we study the relationship between the 
average width (Nu) of dark periods and the detuning 
A,.. We start from investigating the transition rate 
from subspace ^{jOg), to subspace i3{|0e), |le)}. 
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FIG. 7: (Color online) Transition rate 1/Td versus the de- 
tuning Ar. Td represents the average lifetime of the dark 
periods as shown in Fig. 6. The red solid lines show the an- 
alytical results obtained from Eqs. p4j|T^ and the green dots 
show the simulation results. Each green dot is obtained by 
averaging over N = 4000 dark periods. The parameters used 
in both analytical solutions and simulations are: 710 = rie/4, 
= ri,/4, and (a) A = 0, fi™ = lOFie, (b) A = 0, 
fim = SFie, (c) A = IFie, Om = 5Fie . The positions of 
the three peaks are determined by Ar = — V A^ + Q,'^ — A, 
— A and V A^ -f- Q'f^ — A, respectively. It is noticed that the 
heights of the side peaks are equal for A = and unequal for 
A/0. 



i.e., Td = 1/Td, where Td is the lifetime of the dark 
periods. A simple method is to study the time evolution 
of the population of |le) in the qubit-TLS couphng sys- 
tem based on the effective Hamiltonian . As shown in 
Fig. 6 (Here we have deducted the initial-state preparing 
time Tp), the population of |le) is interrupted by ran- 
dom periods of 'darkness'. Quantum jumps from dark 
periods to bright periods indicate that the system jumps 
from subspace A to subspace B. We make an ensemble 
average over the lifetime of the dark periods. As shown 
in Fig. 7, it is interesting that the transition rates reveal 
three peaks as a function of the detuning A^. In addi- 
tion, we find that the positions for the two side peaks are 
f2„i and —Qm, respectively. 

To understand the underlying physics, we rewrite the 
Hamiltonian pT|) in the subspace A and B, respectively. 
In the resonant case A = 0, the Hamiltonian in the basis 



{\A,) ^ {\0g} + \lg})/V2,\Aa) ^ (|0.g) - A = 
(|0e) -I- |le))/\/2, \Bg) = (|0e) - |le))/%/2} reads (see ap- 
pendix) 



Ha- 



r2c/2 fla/2 

I -an/2 -a/2 -a/2 

' f^c/2 -f^c/2 -I- n,n/2 

a/2 -a/2 Ar-nm/2 

(13) 

from which it is clear that the eigenstates in subspace 
A can couple to the eigenstates in subspace B indepen- 
dently with the coupling strength flc/2. Then the phys- 
ical picture is clear as follows. The energy separation 
between the two eigenstates in each subspace is equal to 
the Rabi frequency Sl^. For the detuning = 0, \Ae) 
is resonant with \Be), and \Ag) is resonant with \Bg); 
for A,. = flm, only \Ae) is resonant with \Bg); and for 
A,. = — r^m, only \ Ag) is resonant with \Be)- Note that all 
coupling strengthes are f2c/2, one would intuitively ex- 
pect the height of the central peak to be twice that of a 
side peak, and it confirms from the results of the simula- 
tion. Furthermore, we utilize the Wilcox-Lamb method — 
to obtain an analytical form of the transition rate T]j. In 
the approximation to the first order, the transition rate 
from subspace A to subspace B has a simple form 



02 



ieA.jeB 



(14) 



where 7 — (2710 + Tie)/2, pi is the probability of 
state \i) (here = PAg = 1/2), and Aij are the de- 
tunings between states \i) (|i) € {|Ae), |Ag)}) and \j) 

e {\Be),\Bg)}) in the forms A^e.Be = ^Ag,Bg = 
Ar, AAe,Bg = Ar - l^m, and A^g.Be = A^ + ^m- As 

shown in Fig.7(a-b), our approximate analytical results 
agree with the simulation results considerably. 

Furthermore, Equation (|14p can be straightforwardly 
generalized to the off resonant case A 7^ in the form 



D 



2nt 



E A2. +^2 



(15) 



-sm^ f , 



where Xij represent the coupling coefficients between 
states |i) and Here \i) e {\Ae) = sin||Og) + 

cosf|l5),|Ag) = cosf|0.g) - sinf|l5)}, and |j) G 
{|Be) = sin f|Oe) cos file), = cosf|Oe)- 

sinf |le)} with a = arctan(r2m/A). The coupling co- 
efficients Xij have the forms (see appendix) XAe,Be 

sin f cos |, XAg.Bg = - sin I cos f, XAg,Be 

and XAe,Bg — cos^ f , respectively. The detunings Ay 

have the forms AAg,Bg = AAe,Be = Ar + A, AAg,Be = 

A,.+A+yf}|7+A2, andA^e.Ss = A^+A-v/fir+A2, 
then the positions of the three peaks are given by A^ = 
-^A2 -hO^-A, -A and v^A2 + n'^-A, respectively. 
Conclusively, there are two main differences between the 
resonance case (A = 0) and the off resonance case ( 
A 7^ 0). Firstly, the positions of the three peaks for 
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FIG. 8: (Color online) Average width (No) of dark periods 
versus detuning A,-. The red solid line shows the analytical 
result obtained from Ea. (|14|l -Ea. (|17|l . and the green dot shows 
the simulation result. The parameters used in both analytical 
solution and simulation are: A = 0, 710 ~ rie/4, f2c ~ rie/4, 
fim = lOFie and t,„ = lO/Fie. 



A = are symmetric with respect to = but unsym- 
metric for A 7^ 0. Secondly, in the case A = the heights 
of the side peaks are equal, because all the coupling co- 
efficients between \i) and \j) have the same value 1/2. 
However, in the case A ^ 0, the heights of the side peaks 
are unequal because of the different coupling coefficients 
Ay . As expected, the above analysis is well demonstrated 
in Fig. 7(b-c), where the agreements between numerical 
simulations and analytical results are clear. 

Then we can estimate the average width {Njj) of the 
dark periods in Fig. 5. Supposing the system is in sub- 
space A initially, the probability for the system residing 
in A at the end of a single run with interval is 



P = exp(-r_Dr„). 



(16) 



We have proposed a scheme to observe macroscopic 
quantum jumps in Rabi oscillations of a superconduct- 
ing phase qubit coupled to a TLS. This scheme provides a 
new tool to characterize the TLS quantitatively. Dark pe- 
riods and bright periods in Rabi oscillations indicate the 
TLS residing in states \g) and |e), respectively. There- 
fore, the paramerer To discussed in Sec. V describes the 
transition rate between \g) and |e) of TLS, which is in- 
duced by the coupling between qubit and TLS. To obtain 
the intrinsic lifetime of the TLS, which is known as in- 
duced by coupling to elastic strain field, we just need 
to decouple the TLS with qubit by setting the detuning 
Ar S> rim- Moreover, in recent experiments, the TLSs 
are observed in other kinds of superconducting qubits 
such as charge qubi t^^'^^i^^ and flux qubit^. We be- 
lieve that the method and model in this article can also 
be generalized to such systems as well as other systems 
with similar energy level structure. 

Furthermore, in the process of quantum computing, 
the quantum computer usually evolves into a highly en- 
tangled state. A quantum jump tends to destroy such 
entangled state and make the result of quantum com- 
puting incorrect, which is known as decoherence induced 
by quantum jumps. However, in this article we found 
that quantum jumps can lead to another type of neg- 
ative effect on quantum computing, that is, the errors 
in reading out the quantum state. For example, once 
the system jumps into the dark periods, the qubit's state 
cannot be read out using the normal methods, although 
the qubit can still evolve coherently between state |0) 
and |1). In this case, from the measurement results, we 
always believe the qubit is staying in state |0). That is 
to say, quantum jumps leads to a readout error in the 
qubit. Since the previous literatures emphasized the de- 
coherence induced by the quantum jumps, we hope the 
discussion in this article can bring new insights into our 
understanding of the effects of quantum jumps on quan- 
tum computing. 



Because the probability P for each run is independent, 
the average width of the dark periods can be expressed 
as 



1 



D 



EZi ^"(1 -P) 1 - exp(-rz5r„ 



(17) 



Both the simulated and analytical results for (Nd) are 
shown in Fig. 8 and they agreed very well. The agreement 
between the numerical and analytical results indicated 
the validity of our model. In experiments, the coupling 
strength is typically 20MHz ~ 100MHz,"'2i,30 jg 
typically IMHz^lOOMHz, and other parameters includ- 
ing Fie, Ar and f2m can be controlled flexibly in a large 
regime to fulfill our theoretical discussions by manipulat- 
ing the biased current or microwave amplitude.^'^ There- 
fore, our scheme is feasible within the current technique. 
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VIII. APPENDIX: DERIVATION OF EQ. (fTT|) 
AND EQ. (fT3l) 

In this appendix, we address the derivation of Eqs. 
PT|) and Under the interaction picture and the 
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rotating wave approximation, the Hamiltonian (jlOp be- 
comes 







\ 



( l]„e-^^72 
rJ™e*^72 

17^6^"^'* r2™e-*^72 

a™e*^72 

(18) 

where A = wio — and = lotls ~ ^lo a-re the detun- 
ings. We make a transformation to the rotating frame 
such that the wave function |?/;) in the original Hamilto- 
nian (flS]) can be expressed as 



where 



1^) = uW) 



/ 1 

e^^* 
e'('^+'^-)* 



\ e»(2A+A.)t 

Then the Schrodinger equation 



(19) 



(20) 



•4i 

can now be rewritten as 



where 



u\t)HiU{t)~inu\t)^^^ 

at 

( r2„72 

fJ„/2 A 
f^c A-l-A^ f7,„/2 
\ f^m/2 2A-|-Ar 



(21) 



(22) 



(23) 



This transformed Hamiltonian is exactly Eq. (jlip . To un- 
derstand the relationship between the lifetime of dark pe- 
riods and detunings as discussed in Sec. V, we rewrite 



the Hamiltonian H' in the basis {|^e), l^e), l-Bg)}. 
In this new basis, we find that the Hamiltonian in the 
resonant case A = is actually Eq. fT^ given by 



-ar72 -l^c/2 -f^c/2 
1772 -r272 A,r + an/2 
1772 -r2c/2 A^-f]„72 




(24) 



where 



1 

3 



/ 1 1 
1-100 
0011 

Vo 1-1 



(25) 



For the general off resonance case A 7^ 0, we have 



Ha-b - y^H'v 

( 



UJAg ^Ag,Be^c ^Ag.Bg^c 

^Ae,Be^c ^Ag,Be^c OJBe 
\ ^Ae,Bg^c ^Ag,Bg^c LOBg 



(26) 



w here coAe = (A + 7»2^ + A2)/2, tj ^g = (A - 
^nl^ + A2)/2, c^ Be ^ (A + Vf^2^ + A2)/2 + A + A„ 

C^Sg = (A - ^nl+A^)/2 + A + Ar, XAe.Be = 
Sm2-CO&2-, AAe,Bg — COS , AAg,Be — — Sm , 

Aao Bq = — sin ^ COS f , and 



/ sin f COS f 



1/ 



COS ^ — sin ^ 



(27) 



with a — arctan(il„i/A). 
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